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Vector boson fusion (VBF) processes at the Large Hadron Collider (LHC) provide a unique 
opportunity to search for new physics with electroweak couplings. A feasibility study for the search 
of supersymmetric dark matter in the final state of two VBF jets and large missing transverse energy 
is presented at 14 TeV. Prospects for determining the dark matter relic density are studied for the 
cases of Wino and Bino-Higgsino dark matter. The LHC could probe Wino dark matter with mass 
up to approximately 600 GeV with a luminosity of 1000 fb~^. 



We also show in this Letter that information about 
production cross sections in VBF processes and the dis- 
tribution of pY in the final state can be used to solve for 
the mass and composition of Xi, and hence the DM relic 
density. The cases of pure Wino or Higgsino Xi i as well 
as the case of a mixed Bino-Higgsino x? are studied. 

We note that the production of squarks {q) or gluinos 
{g) through gluon fusion, followed by cascade decay end- 
ing in the production of DM, is the classic setting for 
DM searches in final states with appreciable missing en- 
ergy, multiple jets and leptons. However, determining 
the content of the neutralino and the masses of the su- 
perpartners without any color charges requires specific 
model dependent correlation between masses of colored 
and non-colored superpartners. In very specific settings, 
it is possible to determine the composition of Xi [Hj, as 
well as the mass of light staus or sleptons [12 . In general, 
the combinatoric background poses a major problem for 
such attempts. 

Recently, experiments at the 8- TeV LHC (LHC8) have 
put lower bounds on the masses of the g and q. For com- 
parable masses, the exclusion limits are approximately 
1.5 TeV at 95% CL with 13 fb"^ of integrated luminos- 
ity [TSHIS] . There are also active searches for the lightest 
top squark (t), and exclusion limits in the m^-m^o plane 
have been obtained in certain decay modes p!6 l fT7 ] . 

A direct probe of the electroweak sector using VBF 
processes is complementary to such searches. A variety of 
possibilities exist for the colored sector (compressed spec- 
tra, mildly fine-tuned split scenarios [18^, non-minimal 
supersymmetric extensions, etc.) with varying implica- 
tions for existing and future searches. Experimental con- 
straints (e.g. triggering) significantly affect the ability 
to probe supersymmetric DM in some of the above sce- 
narios, for example those with compressed spectra. The 
important point to note is that a direct probe of the elec- 
troweak sector is largely agnostic about the fate of the col- 
ored sector and provides a direct window to DM physics. 

The strategy pursued in this Letter will be to investi- 
gate direct DM production by VBF processes in events 



Nearly 80% of the matter of the Universe is dark mat- 
ter (DM) [1 . The identity of DM is one of the most 
profound questions at the interface of particle physics 
and cosmology. Weakly interacting massive particles 
(WIMPs) are particularly promising DM candidates that 
can explain the observed relic density and are under in- 
vestigation in a variety of direct and indirect searches. 
Within the context of i?-parity conserving supersymmet- 
ric extensions of the standard model (SM), the WIMP 
DM candidate is the lightest supersymmetric particle 
(LSP), typically the lightest neutralino (xi)^ which is a 
mixture of Bino, Wino, and Higgsino states. 

The DM relic density is typically determined by its an- 
nihilation cross section at the time of thermal freeze-out. 
For supersymmetric WIMP DM, the annihilation cross 
section depends on the mass of x? and its couplings to 
various SM final states, for which a detailed knowledge 
of the composition of x? in gaugino/Higgsino states is re- 
quired. Moreover, other states in the electroweak sector, 
such as sleptons, staus, or charginos can enter the relic 
density calculation. 

It is important to probe the electroweak sector of su- 
persymmetric models directly in order to study their 
DM connection. The main challenge to a direct probe 
of the electroweak sector at the Large Hadron Collider 
(LHC) is the small production cross section of neutrali- 
nos, charginos, and sleptons [2]. 

In this Letter we explore supersymmetric DM pro- 
duced directly at the LHC using vector boson fusion 
(VBF) processes l3J|4j. This appears particularly promis- 
ing since some of the present authors recently showed 
that VBF production is quite effective in probing the 
chargino- neutralino system [5 . It has also been sug- 
gested that VBF processes might be useful in both Higgs 
boson and supersymmetry studies [6UTQ]. VBF produc- 
tion is characterized by the presence of two tagging jets 
with large dijet invariant mass in the forward region in 
opposite hemispheres. As shown in [5], the requirement 
of tagging jets along with missing transverse energy (pr) 
is very efficient in reducing SM background. 



with 2j + ^ in the final state. Such an approach has 
several advantages. The 2j + pY final state configura- 
tion provides a search strategy that is free from trigger 
bias. This is reinforced as the pt thresholds for triggering 
objects are raised by ATLAS and CMS experiments. 

In order to probe DM directly, the following processes 
are investigated: 



PP -^ Xi Xi jj^ 



xtxijj^ xtxijj 



(1) 



The main sources of SM background are: {i) pp -^ 
Zjj -^ i^i^jj and (ii) pp ^ Wjj -^ li^jj- The former 
is an irreducible background with the same topology as 
the signal. The pY comes from the neutrinos. The lat- 
ter arises from events which survive a lepton veto; (iii) 
pp -^ tt+jets: This background may be reduced by veto- 
ing 6-jets, light leptons, r leptons and light-quark/gluon 
jets. 

The search strategy relies on requiring the tagged 
VBF jets, vetoes for 6-jets, light leptons, r leptons and 
light-quark/gluon jets, and requiring large pr in the 
event. Signal and background events are generated with 
MADGRAPH5 ^ . The detector simulation code used here 
isPGS4 go]. 

Distributions of pt(Ji),Pt(J2), Mjij2^ ^^d pr for back- 
ground as well as VBF pair production of DM are studied 
at v^ = 8 TeV and 14 TeV. In the case of pure Wino 
or Higgsino DM, Xi is taken to be outside the exclusion 
limits for ATLAS' disappearing track analysis [21 and 
thus VBF production of xtxti xtxfi ^nd xtxi also 
contribute. The Xi masses chosen for this study are in 
the range 100 GeV to 1 TeV. The colored sector is as- 
sumed to be much heavier. There is no contribution to 
the neutralino production from cascade decays of colored 
particles. 

Events are preselected by requiring ^ > 50 GeV 
and the two leading jets {J1J2) each satisfying pt > 30 
GeV with |A7^(ji, J2)| > 4.2 and rjj^rjj^ < 0. The pre- 
selected events are used to optimize the final selections 
to achieve maximal signal significance {S/^/S~^^). For 
the final selections, the following cuts are employed: (i) 
The tagged jets are required to have pt > 50 GeV and 
^jiJ2 > 1^00 GeV; (ii) Events with loosely identified 
leptons (/ = e^ja^Th) and 6-quark jets are rejected, re- 
ducing the ti and Wjj -^ Ivjj backgrounds by approxi- 
mately 10~^ and 10~^, respectively, while achieving 99% 
efficiency for signal events. The 6-jet tagging efficiency 
used in this study is 70% with a misidentiffcation prob- 
ability of 1.5%, following Ref. [22]. Events with a third 
jet (with Pt > 50 GeV) residing between r]j^ and r]j^ are 
also rejected; {in) The pr cut is optimized for each dif- 
ferent value of the DM mass. For m^o = 100 GeV (1 
TeV), pT > 200 GeV (450 GeV) is chosen, reducing the 
Wjj -^ Ivjj background by approximately 10~^ (10~^). 

The production cross section as a function of m^o after 
requiring |A7^(ji, J2)| > 4.2 is displayed in Fig. [l] The 
left and right panels show the cross sections for LHC8 and 
LHC14, respectively. For the pure Wino and Higgsino 



cases, inclusive XiXi, xtxf, xf xf , and xfxi produc- 
tion cross sections are displayed. The green (solid) curve 
corresponds to the case where Xi is 99% Wino. The 
inclusive production cross section is ~ 40 fb for a 100 
GeV Wino at LHC14, and falls steadily with increasing 
mass. The cross section is approximately 5 — 10 times 
smaller for the pure Higgsino case, represented by the 
green (dashed) curve. As the Higgsino fraction in x? de- 
creases for a given mass, the cross section drops. For 20% 
Higgsino fraction in Xi, the cross section is ~ 10~^ fb for 
m^o = 100 GeV at LHC14. 
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FIG. 1: 
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Figure |2] shows the dijet invariant mass distribution 
Mj^j^ for the tagging jet pair (ji,J2) and main sources 
of background, after the pre-selection cuts and requiring 
Pt > 50 GeV for the tagging jets at LHC14. The dashed 
black curves show the distribution for the case of a pure 
Wino DM, with m^o = 50 and 100 GeV. The dijet in- 
variant mass distribution for W-\- jets, Z-\- jets, and ti-\- 
jets background are also displayed. Clearly, requiring 
^jij2 > 1^00 GeV is effective in rejecting background 
events, resulting in a reduction rate between 10~^ and 
10~^ for the backgrounds of interest. 

Figure [s] shows the p^ distribution for an integrated 
luminosity of 500 fb~^ at LHC14 after all final selections 
except the p^ requirement. There is a significant en- 
hancement of signal events in the high pr region. 

The significance as a function of Xi mass is plotted in 
Fig. [4] for different luminosities at LHC14. The blue, red, 
and black curves correspond to luminosities of 1000, 500, 
and 100 fb~^, respectively. At 1000 fb~^, a significance of 
5(7 can be obtained up to a Wino mass of approximately 
600 GeV. 

Determining the composition of x? foi" a given mass 
is very important in order to understand early universe 
cosmology. For example, if x? has a large Higgsino or 
Wino component, the annihilation cross section is too 
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FIG. 2: Distribution of the dijet invariant mass Mj^j^ normal- 
ized to unity for the tagging jet pair (ji, J2) and main sources 
of background after pre-selection cuts and requiring pr > 50 
GeV for the tagging jets at LHC14. The dashed black curves 
show the distribution for the case where Xi is a nearly pure 
Wino with m^o = 50 and 100 GeV. Inclusive X1X15 xfxt^ 
Xi xf 5 cind Xi xi production is considered. 
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FIG. 3: The pr distributions for Wino DM (50 GeV and 
100 GeV) compared to W-\- jets and Z-\- jets events with 500 
fb~^ integrated luminosity at LHC14. The distributions are 
after all selections except the fJr cut. Inclusive XiXii Xi Xi i 
Xi xf 5 cind Xi xi production is considered. 



large to fit the observed relic density for m^o mass less 
than - 1 TeV for Higgsinos [23] and - 2!5 TeV for 
Winos. On the other hand if Xi is mostly Bino, the 
annihilation cross section is too small. In the first case 
one has under-abundance whereas in the second case one 
has over- abundance of DM. Both problems can be solved 
if the DM is non-thermal [24 (in the case of thermal DM, 
addressing the over abundance problem requires addition 
effects like resonance, coannihilation etc. in the cross 
section, while the under-abundance problem can be ad- 
dressed by having multi-component DM [25 ). If x? is a 
suitable mixture of Bino and Higgsino, the observed DM 
relic density can be satisfied. 

From Figs. [l]and[3J it is clear that varying of the rate 
and the shape of the pr distribution can be used to solve 



FIG. 4: Significance curves for the case where Xi is 99% Wino 
as a function of m^o mass for different luminosities at LHC14. 
The green lines correspond to 3a and 5cr significances. 



for the mass of Xi ^s well as its composition in gaug- 
ino/Higgsino eigenstates. The VBF study described in 
this work was performed over a grid of input points on 
the F — m^o plane (where F is the Wino or Higgsino per- 
centage in Xi)- The pr cut was optimized over the grid, 
and the pr shape and observed rate of data were used 
to extract F and m^o which was then used to determine 
the DM relic density. 

In Fig. [5) the case of 99% Higgsino and 99% Wino were 
chosen, and la contour plots drawn on the relic density- 



''Xi 



plane for 500 fb~^ luminosity at LHC14. The relic 
density was normalized to a benchmark value ^benchmark, 
which is the relic density for m^o = 100 GeV. For the 
Wino case, the relic density can be determined within 
~ 20%, while for the Higgsino case it can be determined 
within ~ 40%. For higher values of tti^o, higher lumi- 
nosities would be required to achieve these results. We 
note we have not evaluated the impact of any degradation 
in pY scale, linearity and resolution due to large pile-up 
events. Our results represent the best case scenario and 
it will be crucial to revisit with the expected performance 
of upgraded ATLAS and CMS detectors. 

In conclusion, this Letter has investigated the direct 
production of supersymmetric DM by VBF processes 
at the LHC. The cases of pure Wino, pure Higgsino, 
and mixed Bino-Higgsino DM have been studied in the 
2j + pY final state at 14 TeV. The presence of the en- 
ergetic VBF jets with large dijet invariant mass as well 
as the large pr due to DM production have been used 
to reduce SM background. It has been shown that broad 
enhancements in the pr and VBF dijet mass distribu- 
tions provide a smoking gun signature for VBF produc- 
tion of supersymmetric DM. By optimizing the pr cut 
for a given m^o , one can simultaneously fit the pY shape 
and observed rate in data to extract the mass and com- 
position of Xi7 and hence solve for the DM relic density. 
At an integrated luminosity of 1000 fb~^, a significance of 
5<j can be obtained up to a Wino mass of approximately 
600 GeV. The relic density can be determined to within 
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FIG. 5: Contour lines in the relic density-m^o plane for 99% 
Wino (blue dashed) and 99% Higgsino (grey dotted) DMs 
expected with 500 fb~^ of luminosity at LHC14. The relic 
density is normalized to its value at m^o = 100 GeV. 



20% (40%) for the case of a pure Wino (Higgsino) for 500 
fb-^ at LHC14, for m^o = 100 GeV. We note that our 
study does not include the effect of large multiple inter- 
actions at high luminosity operations at the LHC. This 
is a very important subject, but outside the scope of the 
present work, because the final performance will depend 
on the planned upgrade of ATLAS and CMS detectors. 
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